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'sj- ' We examine certain classical continuum long wave-length limits of prototype integrable quan- 

•^ I tum spin chains. We define the corresponding construction of classical continuum Lax oper- 

O ! ators. Our discussion starts with the XXX chain, the anisotropic Heisenberg model and their 

generalizations and extends to the generic isotropic and anisotropic gin magnets. Certain 

. , . classical and quantum integrable models emerging from special "dualities" of quantum spin 

\x , 

}_j ■ chains, parametrized by c-number matrices, are also presented. 
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1 Introduction 

Locally interacting discrete integrable spin chains have been the subject of much interest since 
they cropped up in string theory in the study of the AdS/CFT correspondence [1]. Their 
classical, long wavelength limit provides a connection to continuous cr-models describing 
particular dynamics of the string (references on this subject can be found in e.g. [2, 3]). 

This has lead us to tackle here the problem of formulating the classical continuum long 
wavelength limit of the (simpler) quantum integrable closed spin chains in a way that directly 
preserves integrability. Accordingly we will describe the classical Lax-matrix formulation, 
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including the associated classical r-matrix structure, which consistently yields the classical, 
long-wavelength limit, derived for integrable closed quantum spin chain models (see e.g. 

[2, 3]). 

We shall first describe and implement in detail the general Hamiltonian procedure. Then we 
will tackle a number of specific examples, and explicitly compare with already known results 
from alternative derivations. These identifications will establish the validity of our approach. 
We shall in particular consider the paradigmatic example of the long wavelength limit of the 
XXX spin chain, followed by the anisotropic Heisenberg model and the gin classical magnet. 
We finally consider some more complicated cases where the original quantum i?-matrix used 
to build the spin chain by coproduct is "twisted" by a scalar solution of the exchange 
algebra. The corresponding Hamiltonians will be discussed in general, realizing interesting 
formal connections between different classical integrable models. We shall also briefly touch 
upon the inhomogeneous case where the specific twist matrix will be site-dependent. Some 
technical derivations will be exposed in the Appendices. 

Our motivation for this work is to develop a Hamiltonian approach different in its principle 
from the usual Lagrangian formulation of the long wavelength limit, in order to use in cases 
where the latter cannot be applied. In our approach we start from the Hamiltonian inte- 
grability formulation (quantum R-matrix and Lax matrix) guaranteeing a priori Liouville 
integrability of the classical continuous models through a Lax matrix-classical r-matrix for- 
mulation, provided that some consistency checks be made. On all known specific examples it 
will be checked that it yields the same results as the Lagrangian approach. It is indeed a key 
result that the Poisson structure is the same, in all cases when comparison is available, as 
the canonical structure derived from the long wavelength classical Lagrangian. This thereby 
validates the procedure and allows to use it in more general situations where the Lagrangian 
approach may not be used, in particular as a systematic way to build more general types of 
classical continuous integrable models by exploiting the richness of the algebraic approach. 

2 The general procedure 

In this section we outline the general procedure for obtaining a classical Lax formulation 
from the classical limits of the R and monodromy matrices. 



2.1 Classical limit for the i?-niatrix and the monodromy matrix 

A quantum c-number non-dynamical i?-matrix obeys the quantum Yang-Baxter (YB) equa- 
tion [4] 

R\2 Ri3 -R23 = -R23 -R13 -R12 , (2.1) 

where the labels i = 1,2,3 may include dependence on a complex spectral parameter Aj. 
The auxiliary spaces are in this case loop-spaces Vi ® C(Ai), where Vi are (isomorphic) 
finite-dimensional vector spaces. 

Assuming that R admits an expansion ( "semiclassical" ) in positive power series of a param- 
eter (usually denoted h) as 

i?i2 = 1 ® 1 + hru + Oiff) , (2.2) 

the first non-trivial term arising when we substitute this in (2.1) is order of order two and 
yields the classical YB equation 

[ri2, ^13] + [ri2, ^23] + [^13, ^23] = . (2.3) 

This is the canonically known "classical Yang-Baxter equation". It is not in general the 
sufficient associativity condition for a classical linear Poisson bracket, except when r is non- 
dynamical and skew-symmetric (see e.g. [5]). We shall hereafter limit ourselves to such 
situations.^ 

A quantum monodromy matrix T is generically built as a tensor product over "quantum 
spaces" and algebraic product over "auxiliary space" of representations of the YB algebra 
associated to R. Namely, one assumes a collection operators assembled in matrices Lu, 
acting on "quantum" Hilbert spaces labeled by i and encapsulated in a matrix "acting" on 
the auxiliary space Vi. For any quantum space q they obey the quadratic exchange algebra 
[9, 10, 11] 

R12 Liq L2q = L2q Liq R12 , (2.4) 

where operators acting on different quantum spaces commute. The form of the monodromy 
matrix T is then deduced from the co-module structure of the YB algebra 

Ta = Lai La2 ■ ■ ■ LaN (2.5) 



^The dynamical YB equation is related to Drinfel's deformations of quantum groups, whereas the non- 
skew symmetric equation is associated to reflection algebras (see e.g. [6, 7, 8]) and hence to open spin chains 
which we do not consider here. 



and thus naturally obeys the same quadratic exchange algebra (2.4). In particular one can 
pick L = R, the operators now acting on the second auxiliary space identified as "quantum 
space". This way, one builds closed inhomogeneous spin chains with general spins at each 
lattice site (labeled by (i)) belonging to locally chosen representations of some Lie algebra 
(labeled by i). 

We now establish that T has a classical limit by considering in addition the classical counter- 
part of L, labeled by L'^ which then satisfies the quadratic Poisson algebra, emerging directly 
as a semi-classical limit of (2.4), after setting |[A, B] — ?> {A, B}. It reads 

{LUX,), LUX,)} = K,(Ai - A2), LUX,) LUX,)] . (2.6) 

The quantum monodromy matrix has also a classical limit given by (see also [12, 13]) 

^a,{j} = -^al • • • -^aAf • (2-7) 

The exchange algebra for T'^ takes the form 

{T:,Tn = [rat, t:t^]. (2.8) 

This quadratic Poisson structure implies that the traces of powers of the monodromy matrix 
tr{T'^) generate Poisson-commuting quantities identified as classically integrable Hamiltoni- 
ans. In particular, when T^ depends on a spectral parameter, the auxiliary space is a loop 
space V^C{X). Performing the trace over the finite vector space yields a generating function 
tr(T'^(A)) for classically integrable Hamiltonians obtained by series expansion in A. 

2.2 The long wavelength limit 

The usual presentation of the long wavelength limit, such as can be found in [2, 3], is 
a Lagrangian one where the Poisson structure is obtained from the standard derivation 
of canonical variables using a Lagrangian density. Instead, we will present here a purely 
Hamiltonian version of this limit by defining the long wavelength limit of a hierarchy of 
integrable quantum Hamiltonians based on some affine Lie algebra G. We shall define a 
priori the Poisson structure of the classical variables by imposing classical integrability of 
the long wavelength limit of the Hamiltonian through its associated classical Lax matrix. 
We consider a A^-site closed spin chain Hamiltonian H, initially assumed to be governed by 
nearest-neighbour interaction that takes the form 

N 

H = Y,Hu+,. (2.9) 



The classical, long wavelength limit, is obtained by first defining local quantum states as 
linear combinations of the base quantum states, parametrized by a complete set of k con- 
tinuous variables. The number k depends on the choice of G and essentially k = dim(G'). 
These variables, which can be identified as Euler angles in the simplest case of sl{2), become 
the classical dynamical variables once a suitable Poisson structure is imposed. The bras and 
kets are denoted respectively by {n{l,9k)\ and \n{l,9k)), where / denotes the site index and 
6k denote the set of k angular variables. The condition of "closed" spin chain, essentially 
formulated as A^ + / = /, imposes periodicity or quasi-periodicity conditions on the ^fc's. We 
note that we assume that the base quantum states different only by the fact that they are 
defined in distinct sites, hence the frequently used notation below |n/), instead of \n{l,9k)), 
should not be confusing. 

If one considers nearest-neighbor interactions (local) then one defines the classical, but still 
defined in the lattice, Hamiltonian as 

N 

U = ^Ui{t) , Ui{x,t) = {ni\®{ni+i\Hu+i\ni)®\ni^i) . (2.10) 
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For integrable models, we may similarly define the continuum limit of the full set of commut- 
ing Hamiltonians. In these cases the generic Hamiltonians /j(") of the integrable hierarchy 
are obtained directly from the analytic series expansion around some value Aq of the spectral 
parameter of the trace of the monodromy matrix (transfer matrix) as 

oo 

trT(A) = ^(A - Ao)"i/("^ . (2.11) 

ra=l 

By extension, we define in this case the classical Hamiltonians as the expectation value, over 
the N site lattice quantum state, of H^"^' 

7/(")(x,t) = <...(ni|®(n;+i|... //(") ...\ni)®\ni+i)... . (2.12) 

We next define a continuous limit and take simultaneously the thermodynamical limit in 
which A^ — )■ oo. Accordingly, this is achieved by identifying the lattice spacing 5 as being 
of order 1/A^ and subsequently consider only slow-varying spin configurations (the long 
wavelength limit proper) for which 

k ^l{x) , k+i ^l{x + 5) . (2.13) 

In this limit, the finite "site differences" turn into derivatives. 



Given that (2.12) is applied to Hamiltonians of the integrable hierarchy obtained directly 
from the series expansion of the trace of the monodromy matrix, it is immediate that the 
expectation value procedure goes straightforwardly to the full monodromy matrix T (and 
thence to its trace over the auxiliary space which is altogether decoupled from the quantum 
expectation value procedure). Accordingly, we define first a lattice expectation value 

Ta = ...{ni\® {ni+i\... {Lai ^a2 • • • ^aA^) • • • |^z) ® |^/+l) • • • , (2.14) 

which nicely factors out as 

N 

Ta = Yl{ni\Lai\ni) . (2.15) 

1=1 

Assuming now that L admits an expansion in powers of 6 as 

Lai = I + 5lai + 0(5^) , (2.16) 

we consider the product (setting {ni\lai\ni) = la{xi)) 

N oo 

Ta = l[(^ + Sla^ + Y.6-l^:^). (2.17) 

j=l n=2 

Expanding this expression in powers of 6, we get 

i i<j i 

These, multiple in general, infinite series of the products of local terms, are characterized 
by two indices: the overall power n of 6, and the number m of the set of indices i (that 
is the number of distinct summation indices) over which the series is summed. Note that, 
in the T expansion one always has n ^ m. The continuum limit soon to be defined more 
precisely, will entail the limit 5 — ?■ with 0{N) = 0{l/6). We now formulate the following 
power- counting rule, that is terms of the form (for notational convenience lai = lai below) 

m 

^" E ^&^•••t^ E-^=^' (2-19) 

il<i2<---im i=l 

with n > m are omitted in the continuum limit. The latter is defined by 






dx la{x) (2.20) 

and similarly for multiple integrals. Here A is the length of the continuous interval defined 
as the limit of N6. In other words, contributions to the continuum limit may only come 
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from the terms with n = m ioi which the power S"' can be exactly matched by the "scale" 
factor N"^ of the tti- multiple sum over m indices i. In particular, only terms of order one in 
the 6 expansion of local classical matrices Lai = (^j|-^ai|^j) will contribute to the continuum 
limit. Any other contribution acquires a scale factor (5""™ — ). 0, when the continuum limit is 
taken. This argument is of course valid term by term in the double expansion. Being only 
a weak limit argument, it always has to be checked for consistency. 

Let's remark that if L is taken to be R, one naturally identifies 6 with the small parameter 
h, thus identifying in some sense the classical and the continuum limits. However, this is not 
required in general. It is clear to characterize separately both notions in our discussion as 

classical limit : R = 1 + hr , 
continuum limit : L = 1 -\- 61 . (2-21) 

Recalling (2.13), the continuous limit of T, hereafter denoted T, is then immediately iden- 
tified from (2.15), as the path-ordered exponential from a; = to a; = A 

r=Pexp( / dxl{x)) , (2.22) 



where suitable (quasi) periodicity conditions on the continuous variables 6k{x) of the classical 
matrix l{x), acting on the auxiliary space V ^C{X), are assumed. Of course the definition of 
a continuous limit requires that the L-matrices are not too inhomogeneous (e.g. L-matrices 
at neighbor sites should not be too different. This is in fact assured by the long wavelength 
limit assumption. 

2.3 The Lax matrix and r-matrix formulation 

The above identification of T also defines it as the monodromy matrix of the first order 
differential operator d/dx + l{x). In addition, it has been built so that to straightforwardly 
generate the classical continuous limit of the Hamiltonians in (2.12) from the analytic ex- 
pansion 

oo 

tr(r(A)) = ^(A - Ao)"?/(") . (2.23) 

n=l 

We thus characterize l{x) as a local Lax matrix yielding the hierarchy of continuous Hamil- 
tonians 'H^^\ In order for this statement to agree with the key assumption of preservation 
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of integrability we are now lead to require a Poisson structure for I (inducing one for the 
continuous dynamical variables 6k{x)) compatible with the demand of classical integrabil- 
ity of the continuous Hamiltonians. Indeed, such a structure is deduced from (2.6), as the 
ultra-local Poisson bracket 

{/i(x, Ai), hiy, A2)} = [ri2(Ai - A2), h{x, Ai) + hiy, \2)]S{x - y) , (2.24) 

where r is the classical limit (2.2) of the i?-matrix characterizing the exchange algebra of the 
L-operators. More specifically, recalling that Lai = 1 + Slai + 0{6^), plugging it into (2.6) 
and assuming ultra-locality of Poisson brackets one gets 

{la^, lbj} = [rab, Li + hj]^ ■ (2.25) 

One then identifies, in the continuum limit 5 — )■ 0, the factor 6ij/S with S{x—y). Reciprocally, 
it is a well known result (see, for instance L.D. Faddeev's Les Houches Lectures in 1982) 
that if l{x) has a such an ultra-local linear Poisson bracket (2.24) the full monodromy matrix 
between and A of d/dx + l{x) has the quadratic Poisson bracket structure (2.8), thereby 
guaranteeing Poisson commutation of the Hamiltonians. 

We thus obtain a hierarchy of classically integrable, mutually Poisson commuting Hamil- 
tonians from the explicit computation of the monodromy matrix t(A) of the Lax operator 
d/dx+l{x) as /j(") = ^t(A)|A=A()- Such Hamiltonians are however generally highly non-local 
and not necessarily very relevant as physical models. We shall thus extend our discussion to 
local Hamiltonians. 



2.4 The case of local spin chains 



Local spin chain Hamiltonians are more interesting, physically meaningful and easier to 
manipulate. In particular, they are the most relevant objects in connection with string theory 
and the AdS/CFT duality [1]. Their construction generically requires the determination of 
a so-called "regular value" Aq of the spectral parameter such that Lai(Ao) oc Vai, where V 
is the permutation operator. In this sense the expansion of L can be expressed up to an 
appropriate normalization factor as (see also Appendix C) 

L{\) = f\\){l + 51 + 0{5^)) . (2.26) 

Of course only when the auxiliary space a and quantum space i are isomorphic has this 
"regular value" any relevance. One then defines the local Hamiltonians as (denoting as 



usual t(A) = tTaTa{\)) 



^^"^ = ^Mt(A))L^,„, (2.27) 



implying that they are no more such Hamiltonians expressed as linear combinations of higher 
derivatives of t(A). Their long wavelength limit (e.g. (2.10)) is not obviously derivable from 
a straightforward "diagonal" expectation value of the T-matrix contrary to (2.12), since in 
general {F{A)) ^ F{{A)), for any functional of a set of operators A. However, we show 
below that this is indeed the case due to locality properties. Let us first focus for simplicity 
(but, as we shall see, without loss of generality) on the first local Hamiltonian 

i7W = t(Ao)-^^t(A)|,^,^, (2.28) 

where, t~^(Ao) = 'Pi2'P23 • • -Vn-in- This operator acts exactly as a one-site shift on ten- 
sorized states, identifying of course site labels according to the assumed periodicity, i.e. 
A^ + 1 = 1. (Normalization issues will be discussed in Appendix C). Computing the expec- 
tation value of H^^^ we obtain 

(/J(i)) = (nil ® . . . ® iji^\t-\\^) — {^f{\)Tra \{{l + 5la^ + 0{6^))) \n^) ® . . . ® |n^X2.29) 

One has 

(nil ® (nsl ® ... (n^|t~^(Ao) = (risl ® {n^\ ® ... (ni| (2.30) 

and of course A^ + 1 = 1. 

Taking into account the power- counting rule described in section 2.2 we obtain (see also 
Appendix C) that 

TV , N 

(i/«) = J](n,+i|n,)-(r(A)tr,n(l + ^( ^» ) + 0{5^))) . (2.31) 

i=l 1=1 

We then easily establish that in the continuum limit, using the power counting rule and the 
factorized form of both the state vector as (ni| ... (nAr| that the operator to be valued 
over it t'^(Ao) = V12V23 ■ ■ -Vn-in, (^"^(Ao)) = (t(Ao))"^ We finally obtain that in the 
continuum limit 

i^'"') = r'(^o)^t(A)|,^,; = (t(Ao))-^(t(A))|,^,^ = ^(ln(t(A)))|,^,^ . (2.32) 



The computation may be easily generalized along the same lines for any higher Hamiltonian. 
Higher local Hamilltonians are indeed obtained from (2.27), admitting thus an expansion as 

ff(") = t-i(Ao)^t(A)|,^^ + polynomials , (2.33) 

depending only on lower order local Hamiltonians. When computing the expectation value 
of such higher Hamiltonians one gets the expectation value of t~^(Ao)^t(A)|Ao which in the 
continuum classical limit yields 

(^"'(^o):^t(A)U„) = (t(Ao))-^:^(t(A))U„ , (2.34) 



n 



using the same arguments as in the n = 1 case. In addition, one obtains expectation values 
of the polynomials of order k in the local Hamiltonians. In this case expectation values by 
tensor product of local vectors (ni| . . . {tinI are exactly factorized over products of k local 
monomials /ij^ . . -hk,., except if indices i coincide (or at least overlap for multiple indices). 
Locality of the lower Hamiltonians plays here a crucial role. It is clear that such families 
of terms with coinciding or overlapping indices correspond to a second "label" M = k — 1 
and therefore their contribution will necessarily be suppressed in the continuum limit, with 
respect to the contribution of the generic terms (non-coinciding indices) with M = A; by the 
power-counting argument. Hence, it is consistent to conclude that in the continuum limit 

(Polynomial in {H^'^)) = Polynomial in {{H^''>)) (2.35) 

and therefore 

(2.36) 



'dY 



^^"^) = (7T^Mt(A)), ) = :jr,Hm)) 



A=Ao dy- 



A=Ao 



This is the final, key result in systematically establishing the classical continuum limit of 
integrable spin chains. We may now apply this general procedure to all sorts of examples, 
starting with the simpler applications. 

3 The XXX chain 

The XXX model Hamiltonian describing first neighbor spin-spin interactions is given by 

1 ^ 
^ = 2 5Z (^i ^I+i + ^i-i+i + ^I^l+i) • (3-1) 
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It is well known that when one considers the long wavelength limit one obtains a classical 
(T-model [2, 3]. We shall briefly review how this process works. The coherent spin state is 
parametrized by the parameters x, t via the fields 6, if as 

\n{x,t)) = cos9{x,t) e''^^"'*) |+) + sin^(a:,t) e^^'^^"'*) |-) , (3.2) 

where the ranges of variables is ^ G (0, 7r/2) and (/? G (0, tt). One can verify the completeness 
relation 

dfi{n)\n){n\ = 1 , (3.3) 



where the integration measure is given by 

4 
dfi{n) = — sin 9 cos 9 d9 dip . (3.4) 

n 

Then as was described in [2, 3] and in subsection 2.1, one obtains a classical Hamiltonian 
via the expectation value procedure by employing (2.10). The appropriate XXX 2-site 
Hamiltonian is 

Hu+i oc (Vu+i - I) , (3.5) 

where V is the permutation operator acting as V{a ^b) = b^ a ior a, b vectors in V. From 
the definition of "H we are led to compute quantities of the type 

{a\ ® {b\ V \a) ® \b) = {a\b) ® {b\a) = \{a\b)\^ . (3.6) 

They are expressed in terms of scalar products of the form 

(n|n) = cos(^ — 9) cos(v? — (p) + i cos{9 + 9) sm{ip — (p) . (3.7) 

In the long wavelength limit, \n) — \h) = \Sn), 9{x) = 9{x + S) and ip{x) = ip{x + S). We 
conclude that 



H oc f dx (9'^ + sin2(2^) ip'^) . {31 



We shall now derive the Lax representation yielding (3.8) following section 2. The i?-matrix 
for the XXX model is [14] 

R{\) = \ + ihV, (3.9) 
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where for normalization issues we refer to Appendix C This i?-matrix is a solution of the 
quantum YB equation [4]. It has a consistent normalized classical limit defined as 

r(A) = Ip , (3.10) 

which satisfies the classical YB equation. Alternatively, the classical r-matrix may be written 
as 

Set first 

LaniX) = RaniX - j) (3.12) 

and demand that L satisfies the fundamental algebraic relation 

Rab(yXi — A2) Lan{Xi) -Lf,„(A2) = -^611(^2) -^an(Al) Rab{Xi — A2) , (3.13) 

where as usual in the spin chain framework we call n the quantum space and a the auxiliary 
space. Following the general derivation of section 2 and going directly to the continuous limit 
we disregard higher powers in 6 = h (in this case the two small parameters are naturally 
identified). We next define a "local Lax matrix" as a mean value of L on the same coherent 
spin state, taken solely over the quantum space 

{n\Lani\)\n) = l + m{x,\) , (3.14) 

where 

f\{n\a'\n) {n\(T'\n) \ _ 1 f cos2^(x) sin2^(x) e-^^^^^A 

~ \{n\(r+\n) -\{n\(r'\n) J ~ 2 [sin 29 (x) e^^'^^'"'^ -cos29{x) J ' 

where we have used the form of the coherent states to compute the matrix elements explicitly. 
Then / satisfies the classical fundamental algebraic relation 

{h{x, Ai), kiy, A2)} = [ri2(Ai - A2), /i(Ai) + l2{X2mx - y) . (3.16) 

Setting /(x. A) = 11/ A and taking into account the above algebraic relations we get 

{Hi, 02} = Pi2(n2 - ni)5(x - y) . (3.17) 
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The parametrization in terms of the continuum parameters 6{x), (f){x) gives rise to the 
classical version of 5/2- Indeed, parametrizing the generators of the classical current algebra 
as 

S' = cos 29 , S^ = - sin 29 e^^^"^ . (3.18) 

we obtain from the fundamental relation that 

{S+,S-} = S'6{x-y) , {S',S^} = ±2S^6{x-y) . (3.19) 

The continuum parameters 9{x) and (f){x) can also be expressed in terms of canonical vari- 
ables p and q as 

cos26'(a;) = p{x) , ip{x) = q{x) and {q{x), p{y)} = iS{x — y) . (3.20) 

The /-matrix in (3.15) coincides obviously with the potential term in the Lax matrix of the 
classical Heisenberg model. Precisely, one recalls that one must consider as classical Lax 
operator a la Zakharov-Shabat L = d/dx + l{x). The monodromy matrix for L is well 
known now to yield the classical Hamiltonians including the first non trivial one (see [12]) 

(3.21) 

Recalling the expressions (3.18) and substituting in the expression above we obtain the 
Hamiltonian (3.8), hence the process above works consistently. 

Having exemplified the general construction of Section 2 to a simple system and checked the 
consistency of the approach we now turn to more complicated systems by first moving to 
trigonometric and elliptic s/(2) /^-matrices, corresponding to the XXZ and XYZ spin chains. 

4 The anisotropic Heisenberg model 

Consider the generic anisotropic XYZ model with Hamiltonian 

1 ^ 

^ = 2 E (-^-^f ^1+1 + Jy^>'+i + -^^^l^l+i) • (4-1) 

For the following computations it is convenient to set 

Ji = l- 6^a^ , e e {x, y, z} . (4.2) 

13 




The Hamiltonian is written as 

2 ~2 



N 



H-Y.nn.-^-'-^Y. 



i=i 



Af 



O'^a'^i :^ 2, "^^"^ 



z z 



(4.3) 



The additive constant may be omitted here. Taking into account equations (3.5)-(3.8), (4.3) 
and keeping terms of order 6^ we get 

H oc dx fe"^ + sin2(2^) ip"^ + a^ sin2(2^) cos2(2v9) + ay sin2(2^) sin2(2v9) + a, cos2(20)) (4.4) 

This may be seen as a "deformation" of the classical Heisenberg Hamiltonian. The last 
three terms are essentially potential-like terms. In the special case of the XXZ model the 
terms with coupling constant a^, ay are zero, whereas in the XXX case all potential terms 
vanish and one recovers the Hamiltonian (3.8). If we now recall the parametrization (3.18), 
then the expression above reduces to the Hamiltonian of the Landau-Lifshitz model or the 
anisotropic classical magnet [12] 



H (X dx 



dSz 
dx 



dSx 
dx 



dS, 



dx 



"T axO„ -\- ayO -\- azOz 



^y^y 



(4.5) 



We now derive the classical /-matrix for the anisotropic cases. We focus in more detail on 
the XXZ i?-matrix 



R{\) 



sinh(A + f a^ + 



ifj.\ 



smh.{ifi)a 



sinh(?/i)o"+ sinh(A — ^a^ + 



t/iy 



(4.6) 



2 '2 

The classical limit of the XXZ i?-matrix, after appropriate normalization, is given as (we 
divide with the constant factor sinh A) 



i?(A) = l + ifi r(A) + Oif?) 



(4.7) 



where 



r(A) 



1 nf + i) coshA 
sinh A I a+ 



a 



^-^ + i) cosh A 



(4.8) 



The associated classical Lax operator is again obtained from L{X) = i?(A — y) as (once again 
moving immediately to the continuous limit) 



{n\L{\)\n) = 1 + i/i l{x, A) + 0{i/) 



(4.9) 
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where 

/(A) 



1 ( (^1^1^) cosh A {n\a \n) 



sinhA \ (n|o-+|n) -(n|^|n) cosh A 

1 /i^^ cosh A 5- , ^^^^^^ 



sinhA \ 5+ -\S' cosh A 

where S*^, S*^ are the classical generators of the current s/(2) algebra realized in terms of the 
angular variables in (3.18). The continuous variables x,y were omitted here for simplicity 
and will be from now on whenever there is no ambiguity. 

Let us also briefly characterize the classical algebra underlying the model. We set 

, , _ coshA 1 cosh A 1 (a^^\ 

sinhA * sinh(A) * ' sinhA sinh(A) 

Substituting this expressions to (3.16) and taking into account that 

[Ai2, Ai] = -[Di2, A2] , (4.12) 

we end up with the following set of Poisson structures 

{D,, D2} = , {D,, A2} = [Vu, A2]5ix - y) , {Ai, A^} = -[A2, D^]5{x - y) (4.13) 

which give rise to the s/2 Poisson algebra (3.19). 

The full XYZ classical r-matrix also yields, through this process, the classical Lax operator 
of the fully anisotropic classical Heisenberg model, satisfying also the fundamental linear 
algebraic relation (3.16) (see also [12]). A detailed presentation of this derivation is omitted 
here for the sake of brevity. 

5 The gin classical magnet 

In this section we further extend our analysis to the case of higher rank algebras. In partic- 
ular, we study the classical limit of isotropic and anisotropic gin type magnets. 

5.1 The isotropic case 

First consider the generic situation of the isotropic gin quantum spin chain. The i?-matrix 
is given by the general form (3.9), where the permutation operator is of the form 

n 

"P = ^ e-ij <S) Cji , {eij)ki = 5ik5ji . (5.1) 
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The coherent state is parametrized by n continuum parameters as 

n 

\n{x,t)) = '^ai{x,t) Id) , 



(5.2) 



j=i 



where |ej) is the n column vector with one at position i and zero elsewhere. In addition 



(n|n) = 1 ^ ^ |ajp = 1 



(5.3) 



i=l 



Following the process described in the previous sections we end up with the classical r and 
/ operators defined as (in here L{\) = R{X), instead of (3.12)) 



1 1 1 

r(A) = -V , /(A) = T X] <^ij ® {n\eji\n) = T J^ ^^i ^j ■ (5-4) 

The /-matrix satisfies the linear algebraic relation (3.16), which clearly gives rise to the 
classical cuvTent-gln exchange relations among the elements lij{x). These are given by 



{hj{x), kiiy)} = {Siiljk - Sjkhi)S{x - y) 



(5.5) 



We compute next the first local classical integral of motion starting from the spin chain 
Hamiltonian 



TV 






(5.6) 



where we have dropped from the beginning the constant compared to (3.5). Then, defining 
first the Hamiltonian density as 



7/(°)| 



X) 



{n\ (X) (n 



V \n) ® \n) = Y^ hjix) lji{x + 6) . 



(5.7) 



Expanding appropriately this, we conclude that 

1 2 Y^ dlij{x) dlji{x) 



n^'\ 



/ J ''ijK-^) hiy-^' 



i,j=i 



2 ^^-^ dx dx 
«J=1 



(5.^ 



where we have dropped boundary terms by imposing appropriate boundary conditions. The 
first term above is the quadratic Casimir and can be dropped. The second term, proportional 
to 6^, provides, upon integration, the classical Hamiltonian 



(5.9) 
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The classical integrals of motion on the other hand are obtained from the monodromy matrix 
of l{x) for the generic case along the lines described in the Appendix A (see also [12]). 
Comparing with (5.8) they are seen to coincide. 

The direct computation from the classical l{x) matrix is actually presented in the Appendix 
for another model, but it goes along the same lines for the generalized Heisenberg model, 
and is omitted here for brevity. 

5.2 The anisotropic case 

Consider now the anisotropic case. Recall the classical r-matrix associated to Aj^\ [15] 

r(A) = ^^ J2 ^- ® ^n + -^ E e(-"(^->(^-)-^)^e,. e,. . (5.10) 

smh A -^^ smh A ^^-^ 

The associated classical /-matrix will be of the form 

= '-^ Y: ^nisc)e. + -^TTT E e(-"(^-)-(^-)-^)^/,.(a:)e., (5.11) 

smh(A) -^^ smh(A) -^^ 

and satisfies the linear algebraic relation (3.16). Take now the Hamiltonian of the deformed 
spin chain (see e.g. [16]) 

j 
where the matrix f/ is a representation of the Hecke algebra expressed as {q = e^) 

n 

U = E (^*^- ® ""ri - g"''"^'"''^e,, ® e,,) . (5.13) 

It is convenient to rewrite it as 

n 

V = V-l+ ^{l- q'''''^'''^)eu ® Cj, (5.14) 

and also set 

A2 2 

^ = 5a , q--sn{i-j) _ I _ ggn(^ _ ^^^^ ^ ^L . (5.15) 

The Hamiltonian in this case is basically a "deformation" of the isotropic case and again the 
first non-trivial terms arises at order 6'^. We get that 



i,jr=i «^i=i «<i=i 
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It clearly provides a generalization of the Landau-Lifshitz modeP. Note that the last term, 
proportional to a, disappears in the case n = 2, given that in this case ei + 62 = I. A 
generalization starting from the elliptic classical r-matrix can be also obtained, but is omitted 
here for brevity. Similarly, this classical Hamiltonian may be again directly obtained from 
the classical / operator (5.11) as is described e.g. in [12]. 

6 Novel "dualities" of integrable models 

In this section we shall investigate certain integrable "duals" of the XXX spin chain and its 
higher rank generalizations, and we shall derive their classical counterparts. They will be 
based on the coproduct structure applied this time to c-number matrices. 

Given an initial quantum i?-matrix, the c-number YB equation reads 

[i?i2(A), f/if/2] = , (6.1) 

where f/ is a particular scalar n x n matrix. Considering, for instance, R to be the Yangian 
i?-matrix, the latter equation is valid for any n x n matrix. In the case of the XXX matrix 
one may take for example the Pauli matrices to write 

aialR,2{X)44 = R12W ■ (6.2) 

Given the relation (6.1) we may always define a new L-operator L12 = U1L12, which obviously 
satisfies (3.13) as long as L also satisfies it. 

Before we proceed with the classical limit of the models, let us first examine the quantum 
local Hamiltonian arising from the L-matrices. For simplicity we choose A = to be the 
regular value and as usual we define this Hamiltonian as 

i7oc^logt(A)|,^„, (6.3) 

where we now have 

t(A) = tro[Lo^(A) . . . Loi(A)] = tro[f/oi?oiv(A) . . . f/oi?oi(A)] . (6.4) 



^Notc that for q = e*^ one obtains a non-Hermitian Hamiltonian. This is not so surprising given that 
higher rank spin chain as well as higher ^„_]^ afhnc Toda field theories with imaginary coupling are non- 
unitary models. Nevertheless, the relevant physical quantities, such as spectrum excitations, exact S matrices 
etc. have been extensively studied. 



The transfer matrix at A = becomes 

t(0) oc tr^p^V^N . . .U^Vm] 

n = Pl2^23 • • • 'Pn-IN ■ 



f/jv...f/mnf/,_i... UiU^ 



N, 



Taking the derivative of the transfer matrix we find that 



dt{X) 



dX 



(xUn ■■■ f/i+i 



dRii+i{\) 



A=0 



nf/,_i...f/if/jv, 



(6.5) 
(6.6) 

(6.7) 



A=o d\ 

where R = V R and we consider here the XXX i?-matrix. Gathering the information above 
we conclude that the Hamiltonian is again local and reads 



(6.8) 



Note that more general "regularity conditions" of the form Labi^o) = UaVab, will similarly 
guarantee locality of the Hamiltonians derived from (6.3). 

An inhomogeneous generalization of this construction is obviously available by using distinct 
solutions to (6.1) at each site of the chain. One starts from a set of solutions 




«rr«l 



[Ru{X), Ur'U^"] = , 



1,2,...,N . 



The quantum Hamiltonians are derived from the monodromy matrix 

t(A) =tro[Lo^(A)...Loi(A)] =tro[f/f^i?oJv(A)...f/o^'^i?oi(A)] 



(6.9) 



(6.10) 



by applying the co-module structure to generate an inhomogeneous transfer matrix by suc- 
cessive tensoring by Uq Roi{\). Note that the same procedure (albeit with shifts over the 
spectral parameter, or distinct choices of representations of the quantum space i, instead of 
twists by a W c-number constant matrix) is used to get inhomogeneous spin chains in many 
examples. 

Local Hamiltonians are again obtained via (6.3) and following the exact steps of the proof 
of the homogeneous Hamiltonian above. We end up with the generic final expression 



(6.11) 



Such Hamiltonians may be interpreted as describing spin chains in inhomogeneous back- 
grounds modifying the nearest-neighbor couplings in a site-dependant pattern. Possibly 
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interesting types of inhomogeneities include a local defect (one single f/j distinct from 1); 
a domain-like defect (f/j = U for i ^ io and f/j = 1 elsewhere) or a "double-chain" effect 
(f/j = U for even i and 1 for odd i). We shall not discuss the general construction of the 
classical limit for such chain Hamiltonians, postponing it for future studies. 

We now move back to the homogeneous case. A first remark is in order here. It follows from 
(6.6), (6.10) that one expects to have conditions on the choice of f/*-*^ and the parametrization 
of \ni) in order to be able to define classical continuum limits. To illustrate this point let us 
discuss in detail the computation of (t^^(O)). From (6.6) one has 

(t-i(O)) = ®.(n,|(J]f/,)-^ ®, \n,) = J](n,+i|(f/,)-V,) . (6.12) 

i i 

It follows that if Ui and \ni) are such that 

U^\n.^ = \n.^ + 5\vi) + 0(5^) , (6.13) 

where S is the same scaling parameter as for the continuous limit of L and \vi) is some vector, 
which can be chosen to be orthogonal to |nj) without loss of generality. If this condition is 
fulfilled the expectation value then has the following form 

N 

{r\0)) = J](l - {6n,\n,) - 6{n,\v,) + 0{6^) , (6.14) 

where the key identifications hold up to order 5^ in the discrete case. Hence (due to the 
power-counting argument) exactly in the continuous limit 



{rii 



(H U^y ®^ 1^*) = (II^^^I^^I^*))"' • (6.15) 



Hence, the technical derivation of section 2.4 will hold, yielding again (t ^(0)) = (t(0)) ^. 
It is to be expected that similar conditions will arise when considering the derivative term. 

We next examine the most general situation associated to the XXX model. Consider the 
generic 2x2 matrix U and its inverse 

t/=f" ^] , U-^ = ^('^ -^] , D = ad-bc. (6.16) 

\ c dj D \—c a I 

Take also into account the form of U a^ U^^ we conclude that the most general 2-site 
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Hamiltonian after applying the U transformation to XXX is given by the following form 



^a+^{bda'-b^a+ + (fa-) (6.17) 



D 

+— o"" ( — aca^ + a'^a^ — c^a' 

To gain further insight we focus on particular examples. Taking, for instance, the XXX 
chain and setting U = a^, the local Hamiltonian becomes 

N 

H^Y. (^I^I+i + ^'^^1 - ^I^I+i) ' (6.18) 

i=i 

with the characteristic flip of sings in front of the a^ (X> o"^ term. Similarly, for ?7 = cr^, a^ a 
minus sign in front of the a^ ® a^ and a^ ® o"^ terms, respectively, is attached. 

A classical limit can be defined for these modified Lax matrix (recall L{\) = R{\ — y)) set 

also 1/12 = Ui Li2- More precisely, for U = a^ we have, after acting from the left and right 

with the coherent state: 

, ,-,,,, , f A + f cos2^ f sin2^e2^^ \ 

{nL{X n)= ^^ .2 2 _ 6,19 

\— Ysm26'e "'^'^ — A + yCos2^/ 
Now consider the rescaling 9 — )> W, in the small h limit and also appropriately rescale 
A — 7- fi^X. This is precisely a realization of the condition [6.13) on U and \n) discussed above. 
The linear limit of the L operator above becomes after setting iOe^'^'^ = ip, — i^e"^**^ = ip 
(L(A) ~ 1 + hl{X)) 

{n\La{X)\n) = l + thl{X) , [(A) = "2 [ ^ ^^ J • (6.20) 

The /-matrix above is nothing else than the classical NLS Lax operator (see also [17] for lattice 
versions NLS). The new spectral parameter is here defined as A = — p^- hence the critical 
value for A becomes infinite in the continuum limit (in agreement with the computations in 
[12] and Appendix A). 

Consider next U = cr^ and recall the parametrization (3.20). Taking a similar limit we get 
cos 2^ ^hp , sin2^ ^ 1 , e^^''^ -^ 1 ± 2ihq , A ^ ^^A (6.21) 

and keep only lowest order terms. The classical Lax operator takes the following form 

^(A) = 2( ^7^ "^ + ^1 . (6.22) 
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The latter is just the classical Lax operator for the harmonic oscillator. Note that the 
classical /-matrices presented above are of the form: /(A) = XA + B. Taking into account 
that / satisfies the linear algebraic relation (3.16) we end up with the Poisson relations 

{Bi, B2} = V{Ai - A2)6{x - y) , (6.23) 

which lead to the following expected canonical exchange relations 

NLS model : {■?/'(a;), ■?/5(t/)} = 5{x — y) , 

Harmonic oscillator : {^'(a^), viv)} = i^i.^ ^ v) ■ (6.24) 

We next discuss the classical limits of these "dual" quantum Hamiltonians. First, consider 
the quantum Hamiltonian which corresponds to the NLS model (i.e. U = a^) and set 

■H{x) = {a\ ® {b\ (a^ ^a^' + ay^ay -a'^ a') \a) \b) , (6.25) 

where sin 2^ — )■ 2h6, cos 2^ — t- 1. The corresponding Hamiltonian becomes 

H . I ,. OH.) ^ j ,. ,i.m) , (6,26) 

which is just the first integral of motion of the NLS model (see [12] for details on the 
computation of the classical integrals of motion), that is the number of particles. Note that 
higher integrals of motion may be obtained from the higher quantum Hamiltonian (higher 
derivatives of log(t(A))) (see also Appendix B, where higher integrals of motion are also 
computed for NLS). Consider now the situation where U = a^ 

H{x) oc {a\ ® (6| ( - (T^ ® a^ + (T^ ® cr^ + cr^ ® a') \a) ® \b) (6.27) 

and taking into account (6.21), the Hamiltonian density becomes 

n oc p\x) + ^ ^ Hoc fdx (^q\x) + ^^) , (6.28) 

which is simply a classical harmonic oscillator type Hamiltonian, and coincides with the 
Hamiltonian obtained directly from the classical continuum model (see Appendix A). A 
similar Hamiltonian is obtained in the case where f/ = o"^. The only difference is a relative 
minus sign between the p and q terms. The considerations above may be generalized to the 
gin case. Choose for instance 



n-l 










U = N en — e„„ also set 










ki-^l , Inn -> -1 , hn "> "^^i , 


Ini -)■ hlpi , 


kj ^ , 


A ^h^\ , 


) 


i, J T^n . 








(6.29) 
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Notice that above we keep only first order terms, recall also that lij are the generators of the 
classical gin (see section 4). Once again we have implemented condition (6.13). 

From the now standard construction 

(n|L„(A)|n) = 1 + ihla{x, A) , (6.30) 

we conclude that the linear Lax operator in this case takes the form 

n— 1 71—1 

[(A) = -i^ yPi Cm + ipi Cinj - A (^ ^ Cii - e„„j , (6.31) 

j=i 1=1 

which is just the generalized NLS Lax operator (see e.g [18] and references therein). It is 
clear from (3.16) that 

{^i(x), i>jiy)} = 5ij 5{x - y) . (6.32) 

By choosing U given in (6.29) we conclude that the relevant 2-site Hamiltonian is 

^2 = 5Z ^^^ ® ^'fe - X] ^rd ® ez„ - ^ ein ® Crd + Cnn ® e„„ , (6.33) 

acting non-trivially on the sites j, j ' + 1 of the spin chain. Computing now the Hamiltonian 
density in the usual procedure, while recalling the expansion defined in (6.29) one has 

n— 1 

nix) = {a\ ® {b\ H2 |a) ® |6) = . . . oc J2^i(^)^i(^) ■ (6-34) 

1=1 

The later provides the total number of particles of the model (see also [18]) 

N = I dx ^V^i(a:)V5(a;) . (6.35) 

•^ i 

Similar transformations may be found in anisotropic models, however the whole analysis is 
quite subtle in this situation, and will be left for future investigations. 
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A Classical local Hamiltonians 

We compute the classical integrals of motion for the classical harmonic oscillator respectively, 
starting from the associated classical Lax operators. Consider the monodromy matrix 



T{x,y) = Pexpi / dx' l{x',t,\) \ 



(A.l) 



'y 
satisfying the first order differential equation 

^^^^ = l{x)T{x,y). (A.2) 

It may be expressed in the following form [12] 

T(x, y) = il + Wix)) e^(-'^) (1 + Wiy))-' , (A.3) 

where W is anti-diagonal, Z diagonal, and both are expanded at A — )■ oo 

Wix) = J2^^,Zix)= J:^^. (A.4) 

m=0 m,=— 1 

Our purpose is to identify the various W^"^^ and Z'^"^^ and hence the associated integrals of 
motion. 

By substituting the monodromy matrix as in (A.3), and setting /(A) = D + A {D, A being 
the diagonal and anti-diagonal part of the Lax operator) we obtain 

dW 



dx 



+ 2D W + WAW-A = 



^ = D + AW. (A.5) 

dx 



Substituting expressions (A.4) in (A.5) we find 



1^(0) 



a^ 



P^ 



1^(1) = (+^g + ^)ay , (A.6) 

iy(^) = -^(.g + |)V, ... 



and 



dx 

dZ^ 

dx 
dZ^^) 



= -2I, 
, (A.7) 

(?^ + ^)I... , 



dx 4 
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where I is the 2x2 unit-matrix. The first non-trivial integral of motion is obtained essentially 
from the tr J dxZ^^^ 

J«oc /"rfx(g2(a;) + ^!M) . (A.8) 

B Higher Hamiltonians 

We focus on the computation of higher charges in the NLS context starting from the corre- 
sponding quantum model examined in section 6. We shall show that the quantity emerging 
from the quantum higher charge is identical with the higher classical charge, that is the 
momentum. This gives an illustration of the statement in Section 2, eq. 2.36 that the 
second quantum local Hamiltonian derived from the quantum R matrix formulation of the 
spin chain, also becomes in the continuous classical limit the second conserved quantity ob- 
tained from the monodromy matrix derived from the classical Lax matrix l{x). Hence the 
construction is consistent. 

Let us compute the quantum higher charge starting from the quantum NLS Hamiltonian 

j 
Define now the so called boost operator as (see e.g. [19] and references therein) 

B = J2jhn^i- (B-2) 

j 

All higher charges in involution may be obtained via the boost operator B as follows 

So the next charge one obtains via (B.3) is of the form 

if(2) oc J2[h,,+^, /i,+i,+2] . (B.4) 

j 

The three site quantum higher Hamiltonian is then 

/i(2) = [h^I, I^h] , h = a''®a'' + ay ®ay -a'®(r' . (B.5) 

It is now straightforward to show that 

/i(2) = cr^ (g) (or^ (g) or?/ + 0-2/ (g, (j^) - a^ ® (a"" O cr'^ -F cr^ O a"") - a' (^ {a^ ® a"" - a"" ® a^) (B.6) 
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Define now the Hamiltonian density as 

n^^\x) = {n\ (g) {n\ (g) {n\ h^'^^ \n) (g) \n) (g \n) . (B.7) 

Recalling the identifications: sin 2^ — > 2h6, cos 2^ — )■ 1, we have 

{n\(T'-'\n) -^ 2he sin2(^ , (n|a^|n) -^ 2h9 cos2(^ , {n\(r'\n) -^ 1 (B.8) 

and taking into account the first non trivial contribution {0{6)), after expanding the differ- 
ence operators between neighbor sites 6'j+i — 9i^ 9{x + 6) — 6{x) (same for ipi), we conclude 
that ( recall also iOe'^''^ = ^p, -iOe"'^''^ = ^) 

7/(2) = 4^2 ^/ ^ ^(x)iIj'{x) - iIj'{x)iIj{x) . (B.9) 

And indeed the second conserved quantity is the momentum of NLS i.e. 

It is clear that similar computations can be done for higher charges and for other models, 
but these are beyond the intended scope of the present work. We simply focus here on a 
simple example the NLS case to further illustrate the consistency of our approach. 

C Local spin chains: normalization factor 

A delicate normalization issue arises in the considerations of section 2.4. Superficially the 
assumption of "regular limit" of the L matrix Lai{Xo) = Vai clashes with the assumption of 
"semiclassical limit" Lai = l^l + Slai- In fact one is considering two different normalizations 
of the same initial R matrix, yielding respectively the semiclassical L^i matrix and the 
regular Lr matrix. They will differ in the simplest case by an overall c-number factor as 
Lc\ = 5(A — Ao)~ Lr- 

The transfer matrix t yielding °^^' ' | ^_^ , with t^^(Ao) = 'Pi2'P23 • • • Vn-in is obtained from 
application of the co-module structure to Lr- Whenever an overall normalization factor /(A) 
is applied to L the "new" T matrix acquires an overall factor /(A)^ and the Hamiltonian 
H^^^ is shifted by a trivial identity operator A^^ log/(A)|A=Ao- C)ne can then regularize the 
whole construction by suitable substraction of this infinite identity operator whilst keeping 
quantum integrability. 
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The remaining problem in this case is the issue of having to consider the continuous hmit 
of the product Yli (1 $>> 1 + Slu) around the regular value Aq since it would appear from 
the definition of L^i = (A — Xo)^^Lr that the generic term of the infinite product is then 
singular. But in fact one is always dealing with formal series expansion around this value, or 
equivalently integrals on an arbitrarily small contour around but not touching this value. It 
is understood that evaluations of derivatives of the In of the transfer matrix at point A = Aq 
generically mean extraction of the leading term in formal series expansion. The 5 — )■ limit 
is thus always defined, and the continuous limit thus computed will generate, as shown in the 
various examples studied, the correct Hamiltonians, providing a definite check of consistency 
of the procedure. 
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